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Introduction 


Background 

A  significant  attribute  of  laminated  composites  is  their  design  flexibility.  The  layers  or  plies  of  a 
laminate  are,  in  fact,  modular  units  which  can  be  selected  to  provide  distinct  material  properties  and 
fiber  orientations.  It  is  possible,  therefore,  to  tailor  the  properties  of  composites  to  meet  specific 
design  requirements.  Relatively  little  work  has  been  done  to  apply  principles  of  aeroelastic  timing 
to  rotor  blades.  Rotor  blades  must  satisfy  strength,  weight,  fr^uency,  and  lifetime  requirements; 
and  aeroelastic  instabilities  must  lie  outside  the  lUght  envelope,  llie  primary  issue  when  it  comes  to 
design  of  composite  rotor  blades  is  the  overall  suitability  of  tailored  bl^es  for  practical  appli¬ 
cations. 

It  is  well  known  that  elastic  couplings  can  have  a  strong  influence  on  rotor  blade  stability.  Thus, 
any  attempt  to  utilize  the  couplings  of  composite  materials  in  rotor  blade  desi^  must  include  a 
concomitant  investigation  of  the  ii^uence  of  the  coupling  on  stability.  One  possibility  for  utilizing 
tailoring  in  rotor  system  design  is  based  upon  the  creative  use  of  extension-twist  coupling  to 
achieve  enhanced  performance.  This  and  other  uses  have  strong  potential  for  breakthroughs  in 
design  methodology. 

Objectives 

In  this  work  our  objectives  included  (1)  the  development  of  models  and  analysis  methodology  for 
rotor  stability  in  hover,  incorporating  all  relevant  nonclassical  stmctural  effects;  (2)  the  validation 
of  the  models  and  analysis  methodology,  including  convergence;  and  (3)  the  conduct  of  limited 
design  studies  with  the  tools  developed  in  order  to  enhance  understanding  of  the  role  of  elastic 
couplings  on  stability.  We  believe  that  aU  these  objectives  were  met  What  follows  is  a  brief 
summary  description  of  the  results  of  our  work.  Details  may  be  found  in  the  published  papers 
listed  at  the  end  of  the  report 

Unique  Features  of  Work 

New  Models  and  Analysis  Methodology  The  aeroelastic  stability  analysis  of  composite  blades 
should  be  performed  as  a  two-stage  operation  for  the  best  combination  of  computational  efficiency 
and  accuracy.  The  first  step  is  to  determine  blade  cross-sectional  stiffnesses.  This  can  be  done 
either  by  means  of  simplifying  the  analytical  txKxlel  so  that  the  stiffnesses  can  be  determined  in 
closed  form  or  by  applying  powerful  dimensional  reduction  techruques  to  the  three-dimensional 
representation,  which  leads  to  a  two-dimensional  finite  element  solution  for  the  stiffnesses  and  a 
corresponding  set  of  one-dimensional  (global  or  beam)  equations.  The  second  step  is  to  use  the 
sectional  stiffnesses  and  the  corresponding  beam  equations  to  conduct  an  aeroelastic  analysis  of  the 
rotor  blade,  obtaiiting  measures  of  the  global  dynamic  behavior  of  the  blade  such  as  fi^uencies 
and  mode  shapes  wi&out  aerodynamics  or  stal^ty  eigenvalues  with  aerodynamics.  An  optional 
third  step  is  to  substitute  the  results  from  the  global  analysis  back  into  approximate  three- 
dimensional  recovery  relations  to  get  pointwise  asymptotically  correct  approximations  of  displace¬ 
ment,  strain,  and  stress.  Fot  simplicity,  this  thud  step  was  not  undertaken  in  our  work. 

In  the  course  of  our  work  several  sectional  analysis  methods  were  developed,  adopted,  and 
studied.  Analyses  developed  include  one  based  on  potential  energy  and  another  based  on 
complementary  energy  by  Dr.  Rehfield.  Codes  that  were  adopted  and  studied  include  Marco 
Borri’s  code  NABSA,  based  on  Giavotto  et  al.  (1983);  Mark  Nixon’s  code  TAIL  based  on 
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Rehfleld  (1985);  Carlos  Cesnik’s  code  VABS,  based  on  the  work  of  Hodges  et  al.  (1992);  and  a 
code  based  on  ±e  work  of  Badir  (1992).  These  codes  were  used  because  they  have  bwn  shown  to 
be  accurate  and  they  were  available  at  the  time  we  needed  them. 

The  global  beam  equations  of  motion  are  those  of  Hodges  (1990)  and  are  solved  directly  from  the 
mix^  variational  statement  developed  therein.  Our  stability  analysis  was  formulated  from  this 
variational  statement  with  simple  mixed  finite  element  shape  functions.  The  finite  element  shape 
functions  are  the  crudest  possible  based  on  the  variation^  statement  which  is  in  the  “weakest 
possible  form.”  For  example,  all  the  unknown  field  variables  are  represented  as  piecewise  constant 
within  the  element  This  allows  for  exact  element  quadrature  by  inspection.  The  resulting 
coefficient  matrices  are  extremely  sparse,  and  taking  ^vantage  of  that  sparsity  leads  to  a  very 
efficient  computational  scheme  fm*  the  nonlinear  equlilibrium  solution,  which  runs  on  workstations 
in  a  few  minutes. 

The  resulting  aeroelastic  stability  analysis  code,  STAB,  accounts  for  all  possible  material  elastic 
couplings  and  uses  an  aerodynamic  model  similar  to  that  of  GRASP  (developed  at  the 
Aeroflightdynamics  Directorate).  Mathematica  (a  computerized  symbolic  manipulator)  was  used 
to  write  extensive  portions  of  the  code.  STAB  is  C£q)able  of  solving  for  the  equilibrium  deflections 
of  a  rotating  blade  and  linearized  dynamics  of  sm^  motions  about  equilibrium  for  the  hovering 
flight  condition.  Both  of  these  solutions  are  needed  for  stability  analysis. 

Validation  Studies  The  analyses  developed  under  this  grant  were  thoroughly  validated.  The  types 
of  validation  include  linear  and  nonlinear  static  deflection  of  isotropic  and  composite  beams,  free- 
vibration  frequencies  and  modes  for  nonrotating  isotropic  and  composite  ^ams,  and  modal 
damping  for  hingeless  rotors  made  of  isotropic  and  composite  materials.  While  not  exhaustive, 
these  sorts  of  validation  studies  ensure  that  all  structural  effects  are  represented  accurately.  These 
validation  studies  were  done  in  an  incremental  fashion  in  order  to  gain  confidence  in  each  part  of 
the  analysis  and  code.  In  every  case  our  results  are  at  least  as  good  as  published  results,  and  in 
most  cases  they  are  noticeably  better.  In  validation  of  the  large-deformation  characteristics  of 
composite  beams,  for  example,  we  were  able  to  calculate  stiffness  properties  from  “first  principles” 
producing  deflections  which  agreed  with  experiment  much  better  than  existing  theories.  We  also 
validated  our  finite  element  c<^e  against  our  transfer  matrix  code.  Note  that  the  transfer  matrix 
method  is  capable  of  calculating  the  frequency  of  each  mode  with  comparable  accuracy,  whereas 
the  finite  element  results  deteriorate  as  the  number  of  the  mode  increases.  Thus,  if  the  discussion  is 
restricted  to  the  lower-frequency  modes,  the  finite  element  results  can  be  used  with  confidence.  As 
part  of  the  validation,  we  also  determined  the  convergence  properties  of  the  mixed  ffnite  element 
method.  It  was  found  that  4  to  8  elements  were  siffficient  for  static  analysis,  while  16  were 
sufficient  for  stability  analysis.  In  our  study  of  the  modes  and  frequencies  of  nonrotating  beams  we 
found  that  for  three  modes  of  any  one  type  (bending  or  torsion,  say)  to  be  within  2%  requires  32 
elements. 

Concerning  stiffness  and  the  design  process,  analysts  would  generally  prefer  to  work  on  raw 
material  and  geometric  properties  whenever  possible.  For  instance,  summing  ply  thicknesses  and 
presuming  that  each  ply  retains  its  nominal  engineering  properties  makes  life  simple.  We  note  that 
this  simplicity  might  produce  unrealistic  results.  Indeed,  a  trend  described  by  Minguet  and 
Dugundji  was  confirmed  and  studied  further  in  our  work  -  namely,  laminate  thickness  meastu^- 
ments  t^en  in  two  different  stages  of  the  curing  process  produce  different  stiffnesses,  and 
physical  reasoning  leads  to  a  still  different  thickness  and  corresponding  stiffnesses.  The  last 
thickness  is  the  hardest  to  obtain  in  the  design  process,  but  it  leads  to  the  most  accurate  stiffnesses. 
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There  is  very  little  in  the  literature  that  can  aid  in  validation  of  predictive  methods  for  natural 
frequencies  of  vibration  for  anisotropic  beams  of  solid  cross  section.  One  study,  Abarcar  and 
Cunniff  (1972),  does  contain  experimental  results  with  which  results  from  the  present  work  were 
conq)ared.  Our  purely  theoretical  results  for  the  ftee-vibration  characteristics  agre^  with  experinsen- 
tal  data  very  well  -  letter,  in  fact,  than  s(xne  of  the  theoretical  results  published  by  the  experiment¬ 
ers  based  on  idendficadon  of  the  beam  properdes  from  experiment! 

Results  from  STAB  also  agree  very  closely  with  the  U.S.  Army  code  PFLT  for  aeroelasdc  stability 
of  isotropic  blades  in  hover.  We  believe  that  STAB  does  a  better  job  of  calculadng  the  larger 
deflecdons  than  does  PFLT.  Indeed,  some  deficiencies  in  the  PFLT  results  turned  out  to  be  due  to 
its  reliance  on  an  ordering  scheme.  Because  STAB  does  not  use  an  ordering  scheme,  it  does  not 
exhibit  these  problems. 

For  composite  blades,  STAB’s  corieladon  with  Yuan  et  al.  (1992)  is  quite  good,  but  its  correladon 
with  Hong  and  Chopra  (1984, 1985)  is  quite  poor  except  at  zero  ply-angle.  Differences  seem  to  be 
both  quandtadve  and  qualitadve.  The  qualitadve  differences  are  especially  apparent  since  their  cases 
n  and  rv  are  said  to  exhibit  a  flap  divergence  for  Ct/c  =  0.10  in  Hong  and  (Zhopra  (1984, 1985). 
This  divergence,  however,  did  not  appear  when  generating  the  current  results.  In  these  cases  their 
trim  algorithm  diverged,  which  may  have  been  due  to  a  combination  of  modeling  errors  in  the 
work  of  Hong  and  Chopra  plus  their  reliance  on  an  ordering  scheme.  Also,  we  note  that  the 
excellent  agreement  between  the  analysis  of  Yuan  et  al.  (1992)  and  ours,  along  with  the  observa¬ 
tion  of  Fulton  (1992)  that  the  large  variations  in  the  damping  due  to  noodifications  only  affecting  a 
small  fraction  of  the  cross  section,  are  strong  indications  that  the  present  model  is  noore  likely  to  be 
correct  on  this  point 

All  essential  nonclassical  effects  were  included  in  our  model;  indeed,  the  3-D  representation  from 
which  the  model  is  derived  accounts  for  all  possible  deformation  of  the  blades  and  all  possible 
elastic  couplings.  The  approximations  are  in  the  development  of  the  1-D  constitutive  model,  which 
cannot  retain  all  of  the  3-D  information  exactly.  The  nature  of  the  asymptotical  approximation  is 
that  certain  short  wave-length  phenomena,  such  as  the  restrained  warping  effect,  are  not  represent¬ 
ed. 

Modeling  and  Analysis  Development 

A  common  design  for  rotor  blade  sections  is  to  utilize  thin-walled  construction.  This  is  both 
weight-efficient  and  in  accordance  with  the  frequently -used  approaches  to  manufacture  with 
composite  materials.  The  code  TAIL  developed  by  Rob^  Hodges  and  Mark  Nixon,  based  upon 
the  early  work  of  Rehfield  (1985),  applies  to  thin- walled  single  cell  blade  box  configurations.  The 
thin-walled  section,  besides  its  practical  usefulness,  is  a  case  where  closed-form  integral  expres¬ 
sions  for  section  stiffiiess  are  available.  From  a  design  and  manufacturing  point  of  view,  the 
single-cell  configuration  is  likely  to  be  preferred  for  creating  aeroelastically  tailored  configurations 
wi£  elastic  coupling  incorporated,  perhaps  by  layup  over  a  hrxiey  comb  shaped  ewe  that  serves  as 
a  male  tool  for  ease  of  manufacture.  Consequentiy,  special  efforts  have  been  made  to  improve  the 
modeling  of  these  configurations. 

Nearly  all  modeling  accomplished  to  date  has  been  based  upon  displacement  formulations  and 
potential  energy  considerations.  It  is  well  known  that  such  approaches  yield  analysis  results  that 
overestimate  stiffness.  This  is  unconservative  in  stability  studies.  (Consequently,  Dr.  Rehfield 
developed  a  new  modeling  approach  base  upon  complimentary  energy  ccHisiderations  which  tends 
to  underestimate  stiffness.  This  new  theoretical  approach  yielded  two  benefits.  The  first  is  that 
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compliances  or  flexibilities  (the  inverse  of  stiffnesses)  are  determined  directly  in  terms  of  closed- 
form  integrals.  Secondly,  the  out-of-plane  warping  displacement  can  be  estimated  from  closed- 
form  expressions  in  a  manner  consistent  with  the  thm-walled  approximation.  This  formulation  was 
present^  to  ARO  in  an  informal  report  (Rehfield  1991a). 

The  original  work  which  uncovered  a  new  range  of  elastic  coupling  types  (Rehfield  1985)  was  a 
linear,  small-deflection  theory.  This  has  been  extended  to  the  geometrically  nonlinear  range  of 
moderate  rotations.  While  not  as  general  as  Dr.  Hodges’  formulations,  this  approximation  yields 
an  extremely  useful  means  of  evaluating  the  geon^tric  stiffness  effects  of  centrifugal  blade  axial 
loads  (Ml  torsional  stiffness  and  yields  consistent  static  stability  equations  (Rehfield  1991b).  Some 
static  stability  (i.e.,  buckling)  results  will  be  presented  in  a  forthcoming  invited  paper  (Rehfield 
1993). 

With  the  out-of-plane  warping  displacement  known  from  complimentary  energy-based  theory 
above,  it  may  be  used  to  create  a  theory  for  beams  with  initial  twist  which  is,  relatively  speaking, 
simple  and  easy  to  use.  It  was  presented  also  to  ARO  in  an  informal  report  (Rehfield  1991c). 

UmitaticHis  of  time,  resources  and  Dr.  Rehfield’s  major  health  problems  prevented  the  integration 
of  these  works  (Rehfield  1991a,  1991b,  1991c)  into  design  studies.  Furthermore,  considerable 
time  and  effort  must  be  devoted  to  create  sufficient  understanding  to  permit  rational  design  in  the 
presence  of  both  elastic  coupling  and  pretwisL 

Together  with  the  above  modeling  issues  which  are  “analysis”  oriented,  a  preliminary  design 
mcidel  for  cross  sections  has  been  developed  which  aids  in  defining  coupl^  configurations  in  the 
design  process.  We  call  it  the  “Ideal  Tailored  Box  McxleL”  This  is  a  nxxlel  specifically  designed  to 
be  used  in  elastic  tailoring  studies  of  rotor  blades  and  high  aspect  ratio  wings.  It  is  very  simple  and 
serves  as  a  closed-cell  counterpart  to  the  “Ideal  I-Beam  Model”  for  bending  about  a  single  axis.  It 
should  prove  very  useful  for  (1)  teaching,  (2)  optimization  and  parametric  studies,  (3)  facilitating 
physic^  understanding,  (4)  providing  (^alitative  trend  information  and  (5)  isolating  independent 
designer-controlled  mechanisms.  It's  development  was  jointly  supported  by  NASA,  and  it  is  fully 
described  in  (Rehfield,  Chang  and  Zischka  1992).  Also  given  is  an  application  to  a  rotor  blade 
configuration  supplied  to  us  by  Mark  Nixon.  This  section  mcxlel  has  the  added  advantage  that  all 
stiffness  and  compliances  are  found  from  simple,  closed-form  expressions. 

Design  Information 

A  survey  of  companies,  which  included  Bell,  Boeing,  McDonnell  Douglas  and  Sikorsky,  was 
perform^  to  learn  about  blade  structural  design  ^proaches  used  in  industry.  What  has  emerged  is 
a  picture  of  individual  company  approaches  which  depend  upon  company  historical  data  from 
previous  designs,  proprieta^  practices,  iteration,  and  static  and  fatigue  coupon  testing.  It  has  not 
been  possible  to  unify  the  difrerent  approaches.  This  is  made  difficult,  in  par^  by  the  ways  that  the 
individual  companies  organize  their  woik  units  and  by  their  use  of  company  proprietary  software 
that  is  often  po^y  dcKumented  and  has  undergone  evolutionary  mcxiifications  over  time. 

Reliance  upon  historical  data  and  flight  test  data  for  design  insures  that  innovation  will  likely  be 
slow  slow  in  coming.  Parameter  ranges  that  depart  far  from  previous  experience  will  be  resisted. 
In  particular,  balanced  layups  which  exhibit  no  elastic  coupling  are  us^  to  reduce  or  eliminate 
manufacturing-related  warping  of  parts. 

We  believe  that  unbalanced  layups  are  the  most  efficient  fcM*  producing  elastic  coupling.  Relatively 
little  is  known  about  failures  or  damage  processes  in  such  configurations.  Some  preliminary  tests 
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run  at  Northrop  2  suggest  that  a  “scissor-prpe”  damage  mode  can  occur.  In  the  absence  of  more 
information  and  data,  we  have  dealt  with  configurations  which  are  similar  to  those  treated 
previously  in  the  literature.  This  has  the  additional  benefit  of  validating  the  analytical  results. 

Design  Studies 

Two  soft  in-plane  rotor  configurations,  R1  and  R2,  were  used  for  aeroelastic  stability  design 
studies.  Rotor  R2  was  an  extension-twist  coupled  rotor  based  upon  the  rotor  of  Hodges  et  al. 
(1987).  Rotor  Rl,  however,  had  two  laminate  designs,  L2e  ([OVC4])  and  L3e  ([C.C-90®,C.(C- 
90®)2,C1),  both  of  which  were  extension- twist  designs. 

Although  the  laminate  for  rotor  R2  was  fixed,  the  laminates  for  rotor  Rl  had  a  variable  ply  angle. 
The  Rl  laminate  designs,  along  with  their  accompanying  box  beam  dimensions  (which  were 
constrained  to  fit  within  the  selected  airfoil),  were  chosen  because  they  gave  realistic  nondimen- 
sional  rotating  frequencies  which  did  not  significantly  vary  as  the  ply  angle  was  swept  from  zero  to 
ninty  degrees.  In  addition,  the  two  laminate  designs  were  considered  to  fairly  represent  the  entire 
class  of  extension-twist  designs. 

For  rotor  R2,  a  sweep  of  the  thrust  level  was  used  to  investigate  the  lead-lag  damping  variation 
with  hovering  flight  condition.  Results  were  calculated  for  various  stiffness  models.  It  was 
demonstrated  that  a  significant  eiror  appeared  for  this  case  (especially  at  high  thrust  levels)  when 
bending-shear  coupling  was  neglected. 

For  rotor  Rl,  sweeps  of  thrust  level  were  made  for  each  laminate  design  for  various  values  of  the 
ply  angle.  These  studies  demonstrated,  as  suggested  by  Hong  and  Chopra  (1984,  1985),  that  the 
elastic  coupling  available  from  laminate  design  was  able  to  noticeably  affect  the  lead-lag  damping 
level  for  most  thrust  levels.  In  general,  however,  these  laminate  designs  were  not  very  sensitive  to 

either  the  direct  shear  effects  or  the  nonclassical  (parasitic)  coupling  of  bending-shear. 

« 

In  general,  the  accuracy  of  the  composite  predictions  was  found  to  depend  on  the  quality  of  the 
cross-sectional  stiffnesses.  The  two-dimensional  analyses  used,  however,  gave  nearly  identical 
results  for  many  cases  because  of  their  high  quality.  The  performance  of  the  “classical”  stiffnesses 
(which  ignore  shear  deformation  effects),  however,  was  poor  for  several  cases. 

Recommendations 

In  spite  of  the  generality  of  the  work  done  under  this  grant,  there  are  certainly  areas  which  deserve 
additional  worlL  Indeed,  it  should  be  emphasized  that  far  more  design  studies  should  be  done  than 
those  we  were  able  to  complete  once  our  model  was  validated.  Even  so,  the  model  itself  could  be 
made  more  applicable  to  realistic  helicopter  design  by  including  “better”  aerodynanoics,  such  as  the 
Peters-He  model  and  forward  flight  effects,  and  by  also  modeling  the  flexbeam  or  yoke  (a  more 
complete  rotor  system  model).  We  also  strongly  favor  the  creation  of  additional  experimental  data 
sets  for  validation  and  ctHielation  of  composite  blade  analyses. 


2Deo,  R.,  private  communication  to  Dr.  Rehfield,  Northrop  Corporation,  Hawthorne, 
California,  1991. 
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It  is  clear  that  we  need  to  emphasize  “synthesis”  now  that  many  “analysis”  issues  are  resolved  by 
this  wOTk  and  others.  Additional  understanding  of  elastic  couplings  on  dynamics  is  still  needed  as 
are  guidelines  for  design  of  elastically  coupled  rotors. 
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